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Using the collective variables method we revisit the estimates of the Ginzburg temperature for 
the Coulomb-dominated models of ionic fluids. We consider the charge-asymmetric primitive model 



r^ • supplemented by short-range attractive interactions in the vicinity of the gas-liquid critical point. 



For this model, we derive the effective Ising-like i;^^-model Hamiltonian expressed in terms of the 
collective variables describing the fluctuation modes of the total number density. We obtain the 
c/3 ■ explicit expressions for all the Hamiltonian coefficients within the framework of the same approx- 

Cd , imation, namely, the one-loop approximation which produces the mean-field critical parameters. 

g 
I I Based on this Hamiltonian, we consistently calculate the reduced Ginzburg temperature tc for 

C ■ both the pure Coulombic model (a restricted primitive model) and the pure solvophobic model 

■ (a hard-sphere square- well model) as well as for the model parameters ranging between these two 

CS| ' limiting cases. Contrary to the previous theoretical estimates, we obtain the reduced Ginzburg 

>' 

m ' temperature for the pure Coulombic model to be about 20 times smaller than for the pure solvo- 

vn ; 

^^ ■ phobic model. For the full model including the short-range attractive and long-range Coulomb 

Tij- i interactions, we show that tc approaches the value found for the pure Coulombic model when the 

o 

m I strength of the Coulomb interactions becomes large enough. Our results suggest a key role of the 

Coulomb interactions in the crossover behaviour observed experimentally in ionic fluids as well as 
p\ . confirm the Ising-like criticality in the Coulomb-dominated ionic systems. 



I. INTRODUCTION 



For the last two decades, much attention has been focused on the issue of critical and 
phase behaviour of fluids with Coulomb interactions. These studies were motivated by 
controversial experimental results that demonstrated the three types of the critical behaviour 
in ionic fluids: Ising-like critical behaviour, classical mean-field behaviour, and a crossover 
between these two regimes (see Refs [l^M])- In order to interpret the results, the ionic systems 
were classified as either solvophobic or Coulombic. The main attention has been paid to 
the criticality of Coulombic systems in which the phase separation is primarily driven by 
long-range electrostatic interactions. 

It is now generally accepted that the critical behaviour of Coulombic systems belongs to 
the universality class of a three-dimensional Ising model. Earlier experiments that supported 
the expectation of mean-field critical behaviour could not be reproduced in later works (see 

n ri 

Ref. [5|, l6[ and references therein). Precise experiments indicate a crossover from Ising 
to mean-field behaviour characterized by an increase of the non-classical region with the 



polarity increase of the solvent 



,l7|. 



Many theoretical and numerical works on Coulombic systems are based on a restricted 
primitive model (RPM), i.e., an equimolar mixture of equisized charged hard spheres im- 
mersed in a structureless dielectric continuum. It is established that the RPM undergoes 



a gas- liquid- like phase transition at low temperature and low density 8l-lll|. Reliab 



ations 
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timates of the location of the critical point have been obtained using simu 
Simulations also strongly support the Ising critical behaviour of the model 

The major part of theoretical studies of the criticality in the RPM is based on the mean- 
field theories and deals with calculations of the reduced Ginzburg temperature to [la, llTI ]. 
According to the Ginzburg criterion 18|, the mean- field theory is valid only when |t| '^ to, 
where t = {T ~ T^/T^ and T^. are the mean- field reduced temperature and the mean- field 
critical temperature, respectively. Generally, the previous estimates of the reduced Ginzburg 
temperature of the RPM suggest the non-classical region in the Coulombic systems to be of 



the same order, or even larger than in simple fiuids |19l42l|. Thus, excluding the possibility 
of mean-field behavior of the RPM, these results fail to explain the experimentally observed 
reduction of the crossover temperature in Coulombic systems. However, a very small value 
of to for the RPM (about 10'^ times smaller than for a simple fiuid model) was found in 



Ref. 22]. It should be stressed that two conditions are to be taken into account to get a 
reasonable result for the Ginzburg temperature. The fist condition consists in the choice 
of "a reference model" exhibiting typical Ising critical behaviour. The second one is the 
necessity to consider both the reference model and the model studied at the same level of 



approximation. In Refs. 



19 



2l| , the values of the Ginzburg temperature were calibrated by 



a typical model of simple fiuids. However, they were mainly calculated at different, although 
often comparable, levels of approximation. 



On the other hand, the critica 



functional integration methods 23 



i ty 



of Coulombic systems has been studied using the 
27j . In particular, several attempts have been made 
in order to derive the effective Ginzburg-Landau- Wilson (GLW) Hamiltonian of the RPM 



24 



25| . In Ref. |2J], the effective Hamiltonian is obtained in terms of a scalar field conjugate 



to the charge density. However, a non-perturbative renormalization group analysis of such 
Hamiltonian does not allow one to make an unambiguous statement on the nature of the 
critical behaviour of the RPM. An attempt to derive the effective Hamiltonian of the RPM 
in terms of the number density field, strong fiuctuating quantity in the vicinity of the gas- 
liquid critical point, was made using the collective variables (CVs) method in Ref. 25(|. The 
analysis of the Hamiltonian coefficients shows that in spite of the long-range character of 
the Coulomb potential, the effective interactions are of a short-range character and describe 
attraction. The form of the Hamiltonian suggests the Ising-like criticality of the RPM. 
However, the numerical estimations of the relevant coefficients were not presented in this 
work. More recently, an effect of the long-range interactions on the Ginzburg temperature 
has been studied on the basis of the LGW Hamiltonian expressed in terms of the field 
conjugate to the order parameter J26|]. The Hamiltonian coefficients are presented therein 
in the form of an expansion in powers of the ionicity measuring the strength of the Coulomb 
interaction, and the consideration is restricted to the second power. The results have shown 
that an increase in the Coulomb interactions leads to a decrease of the t emp erature region 
of the crossover regime which confirms the experimental observations [29]. A reduction 
of the reduced crossover temperature with an increase of the ionicity is also indicated in 



Ref. 27j, though it is weak compared to the experimental data. The limiting system with 
the pure Coulomb interactions, which is the RPM, lies outside the range of validity of the 
perturbative treatments developed in Refs 26|, 



m 



Recently, non-classical critical exponents have been found for the RPM using the hier- 



archical reference theory [28| . However, an issue of the width of the critical region has not 
been addressed in this work. 

Summarizing, we can state that the nature of the non-classical region in the Coulomb 
dominated systems remains of fundamental interest and presents a real challenge. 

The purpose of the present paper is to derive a microscopic-based effective Hamiltonian 
of the ionic model supplemented by short-range attractive interactions and on its basis to 
consistently calculate the Ginzburg temperature for both the pure Coulombic model and 
the pure solvophobic model as well as for the model parameters ranging between the two 
limiting cases. To this end, we use the CVs-based theory developed for the description 



of phase transitions in ionic systems (see Refs. J30|, |3l|). Following the ideas of Ref. 25| . 
we integrate out the variables connected with the charge-density fluctuations and derive 
the effective Hamiltonian in terms of the variables describing the total number density 
fluctuations. In this paper, we find the explicit expressions for all the coefficients of the 
effective 0^-model Hamiltonian at the same level of approximation, namely, in the one-loop 
approximation corresponding to a one sum over the wave vector. This enables us to get 
consistent estimates for the critical parameters as well as for the Hamiltonian coefficients for 
a whole range of the model parameters. The Ginzburg temperature for the pure Coulombic 
model obtained in this way appears to be about twenty times smaller than for a simple fluid. 
We also study the effect of the interplay of short-range and long-range interactions on the 
Ginzburg temperature. 

The paper is arranged as follows. In Sec. II we give some brief background to the CVs- 
based theory for a charge-asymmetric primitive model with additional short-range attractive 
interactions included. Sec. HI is devoted to the derivation of the effective Hamiltonian of the 
model in the vicinity of the gas-liquid critical point. In Sec. IV we calculate the Ginzburg 
temperature for the hard-sphere square-well model, the RPM as well as for the models 
including both short-range and long-range interactions. We conclude in Sec. V. 

II. BACKGROUND 

We start with a classical two-component system consisting of A^i particles of species 1 
and N2 particles of species 2. The pair interaction potential is assumed to be of the following 



form: 

Uafsir) = ct>^^{r) + 4>%{r) + 0ff (r), (1) 

where (pasi^) is the interaction potential between the two additive hard spheres of diameters 
(Tq, and (7/3. Here, (t)^p{r) is the Coulomb potential: (t)^p{r) = qa(li3(p'~'\^), where 0'^(r) = 
l/(er), e is the dielectric constant. The system consists of both positive and negative ions so 
that the electroneutrality condition is satisfied, i.e., J2a=i laPa = 0, where pa is the number 
density of species a, pa = Na/V, V is the volume of the system. The ions of the species 
a = 1 are characterized by their hard-sphere diameter cxi and by their electrostatic charge 
go and those of species a = 2 are characterized by diameter (T2 and opposite charge —zqq 
(go is an elementary charge and z is the parameter of charge asymmetry). Hereafter we 
consider the case ai = a2 = cr. The potential 0ff (r) describes the short-range attraction. 
We specify 0f^(r) in the form of the square-well (SW) potential of the range A and assume 
^SR^^-j _ ^SR^^^'j _ ^SR^^'j _ (pSR^j.'^^ fjj^g system of hard spheres interacting through the 

SW potential with A = 1.5a can serve as a reasonable model for simple fluids. Such a 
system undergoes a gas-liquid critical point which belongs to the universal class of a three- 
dimensional Ising model. 

Using the CVs method, we can present the functional of the grand partition function of 
the above- described model in the form [30(]: 

^K] = / (dp)(da;) exp f -— ^[0'^-^(A;)pk,7vP-k,Jv + 0^(^)Pk,QP-k,Q] 

+i ^{ujk,NPk,N + uJk,QPk,Q) + In 'E.Hsi'^N - iujN, J^Q " iga^g] I • (2) 

Here, the following notations are introduced. pk,Ar and pk q are the CVs which describe 
fluctuations of the total number density and the charge density, respectively: 

Pk,Ar = Pk,+ + Pk,-5 Pk,Q = Pk,+ — ^Pk,-, 

CV pk,a = Pk a ~iPk a dcscribcs the value of the k-th fluctuation mode of the number density 
of the ath species, the indices c and s denote real and imaginary parts of pk,a; CVs un and 
uq are conjugate to pn and pg, respectively, (dp) and (dw) denote volume elements of the 
CV phase space 

(dp) = n dpo.^n dPk,AdPk,A. (dw) = n "^^o^n ^^^ia^^ia 

A=N,Q k^O A=N,Q k^O 



and the product over k is performed in the upper semi-space (p-k,A = PkA5 '^-KA = '^^a)- 
Coefficients (f)^^{k) and (fP{k) are the Fourier transforms of the corresponding interaction 



potentials. We use the Weeks- Chandler- Andersen (WCA) regularization scheme 32, l33(] for 
the both potentials, ((^{r) and (f)^^{r), inside the hard core. 

'^Hsl^N — ^<^N, ^Q — '^la^q] IS the grand canonical partition function of the hard-sphere 
system with the renormalized chemical potentials 

ZVi + Z/2 _ Z>1 - Z>2 

^N = — T- , t'Q 



where z/^ is determined by 

Ua is the dimensionless chemical potential, z/q, = (5^^—^ In ^a and /Iq, is the chemical potential 
of the ath species; /3 = l/ksT is the reciprocal temperature; A^^ = {2'nma(3~^ /h^y/"^ is the 
inverse de Broglie thermal wavelength. It is worth noting that 0a«(^ = 0) is a finite quantity 
due to the WCA regularization. We introduce 

ic^icAff^k = ic^k.Aff^k - Az/AT, (3) 

where Az/jy = z/^r — z/jy hs with z^Af,HS being the chemical potential of hard spheres. 

In order to develop the perturbation theory, we present InSnsi'^Af.HS — ii^at, t^q — '^<ia^Q\ 



in the form of a cumulant expansion 31 1 



lnSHs[...] = $^^$^ Y. 9Jti^")(A:i,...,A;„K„Q...u;k,„,Q 

n>0 j„>0 ki,...,k„ 

Xa;k,„^j,Ar- • •Wk„,Ar'5ki+...+k„, (4) 

where the prime on a;k,Ar is omitted for the sake of simplicity. In Eq. (jlj), the nth cumulant 
9Jln (/ci, . . . , kn) is a linear combination of the partial cumulants 9?lQ-i...a„(fci, • • • , ^n), the 
superindex i„ indicates the number of variables uj\i,Q in the cumulant expansion. For details 
we refer the reader to Ref. 30|. The expressions for '*Xft( and 9Jt2 are as follows: 

9?tf ^ = (Ar)Hs, m^P = 0, 

^f{k) = G2,ns{k), 'mf\k) = q',z{N)ns 4, 

where (A^)hs and G'2,h5(^) are the average number of particles and the Fourier transform 
of the two-particle connected correlation function of a one-component hard-sphere system. 



It is worth noting that dJll = due to the electroneutrality condition. The recurrence 



formulas for Tin are derived in Ref. [30| (see Eqs. (46) in J30|). 

Taking into account (l3])-(jl]) and replacing pk,Af by pk,Af + ^i (^k we can rewrite Eq. ([2]) 
as follows: 

(3 



'-'Us[j^n,rs,^q]C / (dp)(dc<;) exp f Au^po^N - ^^[(p'^^{k)pk^NP-k. 



N 



~ (— i)" 

+<P^{k)pk,QP-k,Q] + i XlK,iVPk,JV + UJKQPKq] + Yl ^J- Yl Yl ^n"\h, ...,kn) 

k n>2 

+l^ki,Q • • • Wk,„,Q ^k,^_^^,N ■ ■ ■ Wk„,Af^ki+...+k„ j , 



i„>0 ki,...,k„ 



(5) 



where 



Az/AT = Aun — (3 



(N) 



HSTsi? 



C = exp 



Au^iN) 



V 



HS 



(0), 



/3(iV) 



nsJSR 



2V 



(0) 



and InSnsf'^Ar.HS) ^q] = 9^o [^n,us^ ^q\- It is worth noting that the Hamiltonian in Eq. ([5]) 
does not include direct pair interactions of number density fluctuations if (j)^^{k) = 0. 

Since we are interested in the gas-liquid critical point, the small-k expansion of the 
cumulants can be considered. Hereafter we will put 



>(o) 



(0), 



m'^'{k)c:^wi'^'{o) + —m'^ 



>(0) 



(6) 



and approximate cumulants for n > 3 by their values in the long-wavelength limit 

mt\k,,...,K)^mt\o,...). (7) 

III. EFFECTIVE HAMILTONIAN IN THE VICINITY OF THE GAS-LIQUID 
CRITICAL POINT 



In this section, based on Eqs. ([5])-([7]), we derive the effective Hamiltonian of the model ([T]) 
in the vicinity of the gas-liquid critical point. We obtain consistently all the coefficients, 
including the square-gradient term, within the framework of the same approximation. 

First, we integrate over CVs u^^jy and ui^^q in Eq. i^ using the Gaussian density measure 



as the basic one 



25l]. As a result, we get 



^ Z/n 



'='Hs['^Ar,HS, J^Ql 



-I— r 1 -I— r 1 f 

'^ H^^^mTI^^^ / ^^P^ exp (-?/[//„, p^,pq]) 



nmf^ , 



nTlf^ 



(8) 



ki,k2,k3 
X'^ki+k2+k3 



where 

-17 r 1 (0) ^ \^ f (.0) , (2) 

— H[Vcf,pN,PQ\ — — «! P0,N — ^ 2.^ [^2 P]i,NP--k,N- + ^2 Pk,QP-k,C 

k 

^ \^ f W , Q (2) , (3) ^ 

~T\ 2-^ V'^3 P]ii,NPk2,NPk3,N + oa^ Pki,AfPk2,QPk3,Q + «3 Pki,QPk2,QPk3,Q^ 

^ S^ f (0) , a (2) 

li Z^ 1*^4 Pki.Af • • •Pk4,Af + 004 Pki,AfPk2,AfPk3,QPk4,Q 

ki,...,k4 
+4a4 Pki,AfPk2,QPk3,QPk4,Q + «4 Pki,QPk2,QPk3,QPk4,Q ) ^ki+...+k4 + • • • (9) 

and the coefficients a„ (n < 4) are given by 

af^ = -Aun, a'^i^ = 0, (10) 



af ^-^_M-^ 



(0)^2^ 



rf)H ' 9Jt(°) 



2 



,(2) _ 1 / W2) gj^r^gj^r^ o(3^3 



n 

A distinguishing feature of the Hamiltonian (l9l) - ( IT3l) . unhke that for the RPM [25], is the pres- 

(3) 

ence of the coefficients an due to the charge-asymmetry of the model. Using the recurrence 



formulae from Ref. 30|, Eqs. flT0|) -f lT3|) can be reduced to the form given in Appendix A. 



A. Gaussian approximation 

Let us consider the Gaussian approximation, which corresponds to Eq. (Q where only 
the terms with n < 2 are taken into account. In this case, the integration over CVs pk,Af 



and pk,Q in Eq. (|8]) leads to the following expression for the logarithm of the grand partition 
function 






2V 

1 + ^~f{k)mf 



i + ^f\k)mf{k) 



The Legendre transform of In S^ yields the Helmholtz free energy in the one- loop approxi- 
mation. The result is 



/3/rpa 
1 



/3Fi 



RPA 



V 



/5/hs - ^ E^''(^) - ^ E^'(^) + ^/5pV^^(o) 



+ l^rY.^Al + Ppf''S2,ns{k)\+—Y,^ri\l + Pqlzpct>''{k) 



2V 



(14) 

k " ■ k '" 

where /,,, denotes the Helmholtz free energy density. Again, the subscript HS refers to the 
hard-sphere system, S'2,hs(^) is the pair structure factor of a one-component hard-sphere 
system 3J]. The one-loop free energy flT4l) coincides with the free energy in the random 
phase approximation (RPA) of the theory of liquids [34 1 . 

From ( TT^ one can find the chemical potential z/^r in the RPA 

^r = -^,Hs - 1, E ^'^(^) - ^ E ^"(^) + /5p^'^(o) 



Vj. 



+ 



-E- 



/3g2^0^(fc) 



(15) 



J, . + Pqlzpcl>C{k) 
and obtain on its basis the mean-field gas-liquid phase diagram. The gas- liquid critical 



parameters of the charge-asymmetric PM 



i3K( 



0) were calculated in Ref. [lO| using the 



above equations and different regularization schemes for the Coulomb potential inside the 
hard core. 

B. Coefficients of the Landau-Ginzburg- Wilson Hamiltonian 



Our aim is to derive the effective Hamiltonian in terms of CVs pk,Ar related to the order 
parameter associated with the gas-liquid critical point. To this end, we integrate out CVs 
Pk,Q in Eq- dH]) following the programme outlined for the RPM 25|. As a result, we arrive 
at the effective LGW Hamiltonian of the following form: 



- W [Pn] - - ^ ,/,rv„_i 2^ Ctn Pki.Af Pk2,Af • • • Pk„,Ar 5ki+...+k„ 
n>l -^^^Z 



(16) 



•i^n 



Coefficient a„ can be presented as 

a„ = a^n^ + Aa„, (17) 

where the second addend is the correction obtained after integration over CVs pk,Q- Each 
correction Aa„ has the form of infinite series. In particular, after some algebra we find the 
following expressions for Aa„ in the approximation corresponding to a one sum over the 
wave vector q 

A«i = -J^W^ E ^w(^) + • • • ' (18) 

q 

Aa2(A:) = ^^J]GQQ(g)-^^l^J]GQQ(g)GQQ(|q + k|) + ..., (19) 

q q 

Aa3(fci,/c2) = -^^^<5QQ(g) + ^^J]GQQ(g)GQQ(|q + ki|) 

q q 

^E^QQ(^)^w(|q + kil)5QQ(|q-k2|) + ..., (20) 



<N > 

q 



AaA{ki,k2,h) = —^T—^GQQ{q) - ——— ^GQQ{q)GQQ{\q + k^ + k2\) 



<N>^^ ^^^^' <N> 

q q 



^^ $^GQQ(g)GQQ(|q + ki|) + -^5^GQQ(g)GQQ(|q + ki|)GQQ(|q-k2-k3|) 



<N>^^ ^^'^' ^^"^ " <N> 

q q 



-^W^ E GQQ(g)GQQ(|q + ki|)GQQ(|q - k^DGggdq - k2 - ksl) + . . . , (21) 

q 

where we have made use of the formulas from Appendix A. Hereafter, for notation simplicity, 
the subscript HS is omitted. 

Ggqlq) is the Fourier transform of a charge-charge connected correlation function of the 
charge-asymmetric PM determined in the Gaussian approximation 

1 



Ggqiq) 



1 + Ppqiz(pc{q) 



For the WCA regularization 33 1 



~^, , , sm OCT .. 

<P^{q) = 47r— ^. (22) 



eq-" 
10 



Now, we are in a position to study Eqs. ( !T8|) -( l2T]) in detail. As is seen from Eqs (fT8l) - (l2T!) . 
corrections Aa„ for n>2 depend on the wave vectors. Expanding the second term of Aa2 at 
small k one can readily see that the linear term vanishes. As a result, Aa2 can be presented 
as follows: 

Aa2 = Aa2,o + A;2Aa2,2 + C(A;^), (23) 

where 

Aa2,o = T^^GQQ{q)-——^\GQQ{q)f, 

^«2,2 = - ^ E GQQ{q)&S'Q^q) (24) 

and the superscript (2) refers to the second-order derivative with respect to the wave vector 
k: GQQ{q) = d'^GQQ{\q^ + 'k\)/dk'^\k=o- The appearance of the factor k'^ in Eq. ( 12^ is caused 
by the charge-charge correlations being taken into account. It should be noted that the 



expression for Aa2,2 has the same form as that obtained for the RPM in Ref. 25|. However, 
the numerical estimates of the effective Hamiltonian coefficients were not made in that work. 
It will be done below. 



The integrals J GqQ^q) dq entering Eqs. (iT8l)- (!2Tl) are divergent at g — )■ oo. This 
divergence can be avoided by introducing the cutoff wave- vector which, however, leads to the 
cutoff-dependent results. Another way is to rearrange the terms in Eqs. flT8l) - fl2Tl) expressed 
by infinite series. Here, we will follow the second way. 

We approximate Aa2 by Eqs. fl23|) - flM|) and replace Aa3(fci, . . .) and Aa4(A;i, . . .) by their 
values in the long- wavelength limit putting Aa3(A;i, . . .) ~ Aas^o and Aa4(A;i, . . .) ~ Aa4,o. 
It is convenient to calculate Aa„ using Feynman diagram presentation. Limiting the series 
in Eqs. (fT8l) -( l2Tl) to the order of one sum over q one can present the part of the effective 
Hamiltonian (TT6|) in the following diagrammatic form 

Aaipo,Ar = ----Q, (25) 

1 1 

2!(iV) ^ ^^-'OP^^^P'KN = D " 2 - O "' ^^^^ 

k 

I 

J— Y^ Aa3,oPki,ivPk2,JvP-ki-k2,JV = ~^ "!^0 + ^ " \3^s ~ y\ ' (27) 



3!(iV) 



11 



JTT^ 2_^ A«4,0Pki,AfPk2,AfPk3,AfP-ki-k2-k3,Af — 12;/ J — 6 ( \ 



4\{N) 

^ ' ki,k2,k3 



(28) 



In Eqs. f l25|) -f l28|) . each external leg represents CV Pk.iv, the vertex represents the factor 

— j-'^ki+...+k„ and each internal line corresponds to the propagator GqQ^q). Integrals over 
wave vectors are implied here. As is seen, the restriction to a one sum over q in the 
expressions for Aa„ leads to the one-loop diagrams having the same number of vertices and 
bonds. 



Now, introducing 



we substitute 



^(g) _ /3pgo#^(?) 



l + /3pg2#^(g) 



GQQ{q) = l+^{q) (29) 

into Eqs. ( IT8l) -( l2T1) . Then, adding all the terms corresponding to the one-loop diagrams in 
Eqs (I25l) -( l28l) we present Aon in the form: 

A«i = -^-7V75Z5^(^)' (30) 

q 

A«2(A;) = -^^^ E \3{<l)f - ^ E 5^^'^(^) [1 + ~9{<i)] , (31) 

q ^ ' q 



A«3 = -^^^E[^(^)]'' (32) 

q 

A«4 = -^^E[^(^)]'- (33) 

q 

The integrals entering Eqs. (l30!) - (l33l) are convergent in the lower and upper limits. Again, the 
superscript (2) denotes the second-order derivative with respect to k. Comparing Eqs. fl3T]) 
and fl2^ it is easy to check that substitution (!29l) does not change the coefficient Aa2,2- 
Finally, inserting Eqs. ( I30|) - (l33l) in ( IT71) we arrive at explicit expressions for the coefficients 



12 



Qn in the one-loop approximation: 

aifl = -Aun - Ci,c, (34) 

«2,0 = -P C2,HS + f^P^^'^iO) - p C2,c, (35) 

a2,2 = -ipCSs + ^/3P0'^'^'^ - 4(^5^5^^'^ (^) [1 +5^(^)] ' (36) 

q 

^3,0 = -P^Cg^HS - P^Cg^c, (37) 

«4,0 = ~P C4^HS ~ P ^4^0) (38) 

where a„_o = «n|fci=0; c^2,o and 02,2 are the coefficients of the expansion 02 = 02,0 + ^^^2,2, and 
C2,Hs(^) = C2,HS + ^^^^2 Hs- ^n,HS is the Fourier transform of the n-particle direct correlation 
function of a one-component hard-sphere system in the long-wavelength limit 

-pC2,Hs = ^^, p'^-^C„,Hs = p"-'C„.Hs(0,...) = «i°\ n>3 (39) 

and Qhs — d'^C2,Hsik) /dk'^\k=o- Using the formulas from Appendix A, one can establish a 
link between the direct correlation functions C„,hs and the connected correlation functions 
G„,Hs for n < 4. 

C„,c denotes the contribution to the ra-particle direct correlation function from the charge 
subsystem 



p"-^C„,c = p"-^C.,c(fc = 0) = ^^^ E l9i<l)r ■ (40) 



q 

It is easy to check that the coefficients a„,o for n > 2 can be obtained from the one-loop free 

energy (see Eq. ([E])) 

_ n-i <9"(-/3/RPA) _ „_i^ . . 

O-nfi — P -K-;i — P (-nlU, ...J, 

where C„(0, . . .) denotes the Fourier transform of the n-particle direct correlation function of 
the full system in the long- wavelength limit. It is worth noting that the functions C„ differ 
from the ordinary direct correlation functions c„ by an ideal term 3J, |35|. 

Coefficient ai,o is the excess part of the chemical potential u^ connected with the short- 
range attractive and long-range Coulomb interactions. Equation ai,o = leads to the 
expression for u^ in the RPA (see Eq. (IT5l)). Coefficient 02,2 describes the square-gradient 
term. We emphasize that all the coefficients given by Eqs (!M|) - (!38|) are found in the one- loop 
approximation corresponding to a one sum over the wave vector. The charge asymmetry 
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does not manifest itself in the above equations and, therefore, there is no difference between 
the RPM and the charge-asymmetric PM at this level of approximation. The resulting 
effective Hamiltonian reads as 

■^ ' k '^ ' ki,k2 

XPki,ArPk2,AfP-ki-k2,Af + a]/ t\t\3 / . '^4,oPki,AfPk2,AfPk3,AfP-ki-k2-k3,Af • 

■^ ' ki,k2,k3 

Therefore, within the framework of the same approximation, we have derived the 
microscopic-based expressions for the coefficients of the effective Hamiltonian. Taking into 
account the charge-charge correlations through integration over the charge subsystem (the 
CVs pk.g) we get a contribution to the coefficients at the second order which describes the 
effective attraction of short-range character. The resulting Hamiltonian has the structure of 
the LGW Hamiltonian of an Ising model in an external magnetic field. Eqs. flM|) -f HQ|) will 
be used below for the calculation of the Ginzburg temperature. 

IV. GINZBURG TEMPERATURE 

We present 02,0 in the form: 

_T-Tc 

0'2fl — ^2,0 + 0,2,t't, t — — , 

J- c 

da2n 



where 02,0 = a2,o(^ = 0) and 02,* 
value. 



. Hereafter, the subscript c refers to the critical 

i=0 



At the critical point, the system of equations 

^2,0 = 0, 03,0 = (41) 

holds yielding the critical temperature and the critical density. After the substitution of T^ 
and pc in the equation ai^o = one gets the critical value of the chemical potential u^. 



Following [l8|, ll9| , the reduced Ginzburg temperature can be written as follows: 



to = :^2^^ (42) 

where all the coefficients should be calculated at the critical temperature and the critical 
density determined from Eqs. ( HTj) . 
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A. Hard-core square- well model 

First, we consider a one- component system of hard spheres interacting through the SW 
potential of depth e and range A = l.So". This system exhibits a typical Ising critical 
behaviour. For the WCA regularization, the Fourier transform of the SW potential has the 

n 

form [32] 

4)^^{k) = 0^^(O)— ^[-Ax cos(Aa;) + sin(Ax)], (43) 

where x = ka and 0"^^(O) = — ecr^^A^. 

For this model Aa„ = 0, and we get from Eqs. flMl) - fl38l) and Eq. fH3l) simple expressions 
for the coefficients of the LGW Hamiltonian 

ai,o = -Au^ + — — - ---, (44) 



22^SR 2^SR 

~ 8r]X^ _ 8r/A3 _ 1 ~(2) 4 r/A^ 

«2,0 - -PL.2,HS - 77^55-5 a2,t-T^;^^, '^2,2 - --PC2.HS + 777;qR' \^^) 



^3,0 = -P Cs^HS, ^4,0 = -P C4,HS, (46) 

where Aun = i^n — ^n,}is- The reduced temperature in the hard-core square- well (HCSW) 
model, T^^, is defined as the ratio between the thermal energy and the interaction strength 
of the two hard spheres at contact T^^ = k-oT /e, and t] = iipa^/G is the packing fraction. 



We use the Percus-Yevick (PY^ 



of the hard-sphere system 3^, l36|. As a result, we have 



approximation for the n-particle direct correlation functions 



~ (1 + 27^)^ ~^2) _rKl6-m7 + V 

-PC2,HS - — ITT' -P^,HS - inn _^^4 



(l-r^)^' ^-^^'^^ 10(1 -r^)^ 

(l-7r/-6r7^)(l + 2r7) 
-P C3,HS - ^3^^y^ , 

3~ _ 2{l - 6ri + 15r]^ + 56ri^ + 2Ari'^) 

-p ^4-Hs jY^:^^ 

Using Eqs. (14T1) . (l42l) and ( l44l) -( l46ll one can easily calculate the reduced Ginzburg tempera- 
ture. The corresponding results are presented in Table 1 (the first row). It should be noted 
that we get tg = 1.55 if we pass to the expression for the the reduced Ginzburg temperature 
used in Ref. |2(| which differs from Eq. (H2l) by the factor (Pc)^^- This agrees with the result 



to = 1-57 obtained in 20|] for the same system in the RPA. 
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B. Restricted primitive model 

Now, we consider the RPM setting 0^^ = in Eqs. ([3lD-([38D. Using Eq. (El]), we get 
the following explicit expressions for the coefficients of the effective Hamiltonian 

ai,o = -Az/iv - ^^ + iu (47) 

^^2-0 = (^_^)4 + ^2, «2,t = ^ I2,t, (48) 

_ 7^(16 -llr/ + V) 1 

"^■^ - 20(1 -r^)^ 6^^'^^' ^^^^ 

(l-7r7-6r/2)(l + 2r^) . 
03,0 = (^_^)5 + ^3' (50) 

2(1 -6ri + 15r/2 + 56r/3 + 2 V) 
04,0 = (l^^^ ^ ''■ ^ ^ 

In Eqs. ( 147|) -( 15T|) the following notation are introduced: 



(n-l)!(-«:2)_ 

'0 



ttT^ 



2^n-l /"oo 

x2 



sm uc 



x-^ + fi;2sin(a;) 



dx, 



°° x^sin2(x) 



^2,t = / jdx, 

Jo (x^ + fi;2gin(x)) 

ii2 = / x^ [k^x^(1 + cos^(x)) + 2(a;'^ — 2^2 sin(x))(2xcos(x) — 3sin(x)) 
Jo 

+x^ sin(x)] /{x^ + K^ sin(x))'^dx, 
where k = Koa with k|, = AirpPq'^z being the Debye number. 

TABLE I: The reduced gas-liquid critical parameters, the coefficients of the effective Hamiltonian 
and the reduced Ginzburg temperature tc for the hard-core square-well model (HCSW model) and 
for the restricted primitive model (RPM) in the one-loop approximation. The superscripts SR and 
C refer to the HCSW model and RPM, respectively. 

Model Tc '^ p* a2,t 02,2 04,0 *G 

HCSW 1.2667 0.2457 2.7426 0.4536 2.7421 0.0937 

RPM 0.08446 0.0088 1.0758 0.2570 0.1752 0.0053 

The reduced temperature T*-^is defined in the standard way as the ratio between the 
thermal energy and the interaction strength of the opposite charged hard spheres at contact 

T^ = ^-^. 52 
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The loci of equations a2,o = and 03,0 = are shown in Fig. 1. The two hnes intersect at 
a maximum of the gas-hquid spinodal yielding the coordinates of the critical point in the 
one- loop approximation (see Table 1). 




0.004 0.008 0.012 0.016 0.020 



FIG. 1: Restricted primitive model: the loci of equations 02,0 = (solid line) and a^fi = (dashed 
line). Temperature T is given by Eq. ()52p and rj = irpa^/G is the packing fraction. 



Substituting the critical temperature and the critical density into Eqs (H8|) - (!5T]) we cal- 
culate the coefficients of the effective Hamiltonian. Then, we get from Eq. ( |42|) the reduced 
Ginzburg temperature of the RPM. The results are presented in Table 1 (the second row). 
As is seen, the reduced Ginzburg temperature found for the RPM is about twenty times 



smal 



mgs 



er than that obtained for the HCSW model. Therefore, contrary to the previous find- 
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2l[ |. our results suggest that the critical region for the RPM is much narrower than 



the critical region for a non-ionic model. 



C. Long-range interactions versus short-range interactions 



Here, we consider the full system given by Eq. ([T]) and study the effect of the interplay 
of short- and long-range interactions on the Ginzburg temperature. We briefly call the 
model as a RPM-SW model. For the RPM-SW model, we get explicit expressions for the 
coefficients of the effective Hamiltonian combining Eqs. ( H^ -( l46l) and Eqs. ( H71) -(l5T]). We 
also introduce the parameter 

a 



YSR 



Qoz 



yc 



eae 
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measuring the strength of the Coulomb interaction with respect to the solvophobic interac- 
tion. 



0.10 




0.020 



FIG. 2: RPM-SW model: the loci of equations 02 = 0, (solid lines) and 03^0 = (dashed line) for 
a = 25, 50, 100 and 00 (RPM). T^ is given by Eq. (f52|) . rj = iipa^ /Q is the packing fraction and 

a = qQz/{e(j£). 



2.5 


/-T— 


0.5 


2.0 


- / \ 


^\. 


1.5 


- / 


1 


1.0 


- / /-""^T^^ 




/ / : 2 


^— — _ 


0.5 






// _--—'''''''" / iU 


n n 


1 . 1 . 1 . 1 . 1 . 



0.05 0.10 0.15 0.20 0.25 0.30 



FIG. 3: RPM-SW model: the loci of equations 02,0 = (solid lines) and 03^0 = (dashed line) for 
a = 0.5, 1, 2, 5 and 10. The meaning of the symbols is the same as in Fig. 2. 



As before, we first solve the equations for the critical parameters. The loci of these equa- 
tions are shown in Fig. 2 for a = 25, 50, 100 and oo (RPM) and in Fig. 3 for a = 0.5, 1, 2, 5 
and 10. As for the RPM, the curve a^Q = intersects the gas-liquid spinodals at a maxi- 
mum. The results for the critical parameters, the Hamiltonian coefficients and the reduced 
Ginzburg temperature are presented in Table 2. As is seen, the reduced Ginzburg temper- 
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ature first decreases with an increase of a and tlien begins to increase slowly approaching 
the RPM value for a > 100. Quite probably, the RPM-SW model with a certain value of 
a lying between 10 and 25 undergoes a tricritical point (see, e.g., Ref. [37|, l38|). This issue 
deserves a separate study. All the coefficients have the same trend with an increase of a 
from to oo, and their trend coincides with the trend of Iq. Both T'^ and p* (or r^c) decrease 
when a increases and reach the critical values of the RPM. The trend of T^^ is opposite to 
the trend of T^ . 



In order to establish a link to our previous results (see Ref. [26|) we pass from a to the 
ionicity X by means of the relation 

a = T^^X, where ^=^- (53) 

Fig. 4 presents the dependence of the reduced Ginzburg temperature on the ionicity. The 
;he present work are shown by solid circles while the line represents the results 



results of 

from Ref. 26|. As is seen, the both groups of results agree well for X < 6.25 corresponding to 
a < 10. The deviations between the results that appear for large values of X are due to the 
different approximations used in calculating the Hamiltonian coefficients in the present and 
previous works. Therefore, the results of this paper generally agree with the previous results 
obtained within the framework of a different perturbative treatment and, in turn, confirm 
a key role of the Coulomb interactions in the reduction of the crossover region compared to 
the solvophobic interactions. 

V. CONCLUSIONS 

In this paper, we have revisited the issue of the criticality of the Coulomb dominated sys- 
tems. The model considered, besides Coulomb interactions, includes short-range attractive 
interactions. Using the CVs-based theory we have derived the 0^-model LGW Hamiltonian 
in terms of the CVs describing the fluctuation modes of the total number density. The result- 
ing form of the effective Hamiltonian conflrms the fact that the criticality of the Coulombic 
models belongs to the universal class of a three dimensional Ising model. The important 
feature of the developed approach is that it enables us to obtain all the coefficients, including 
the square-gradient term, within the framework of the same approximation. 

The LGW Hamiltonian obtained has been used for the calculation of the gas-liquid crit- 
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TABLE II: The reduced gas-liquid critical parameters, the coefficients of the effective Hamiltonian 
and the reduced Ginzburg temperature tc for the RPM-SW model (see the text) in the one-loop 
approximation . 



a 



t9 



rpSK 

c 



Pi 



a2,t 



«2,2 



«4,0 



to 



0.5 


2.56537 


1.28269 


0.2429 


2.7132 


0.4506 


2.6606 


0.0903 


1 


1.29904 


1.29904 


0.2393 


2.6694 


0.4427 


2.5441 


0.0885 


2 


0.66471 


1.32943 


0.2328 


2.5908 


0.4288 


2.3275 


0.0840 


5 


0.28387 


1.41935 


0.2170 


2.4320 


0.3986 


1.8133 


0.0676 


10 


0.16012 


1.6012 


0.2026 


2.4104 


0.3740 


1.2586 


0.0390 


25 


0.09054 


2.2636 


0.0159 


1.0406 


0.2175 


0.0906 


0.0024 


50 


0.08705 


4.3524 


0.0109 


1.0597 


0.2396 


0.1472 


0.0047 


100 


0.08568 


8.568 


0.0097 


1.0676 


0.2481 


0.1602 


0.0050 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 




FIG. 4: The reduced Ginzburg temperature as a function of the ionicity (see the text" 
indicate the results of the present work and the line represents the results from Ref. 
by Eq. dMD. 



'he symbols 



26l | . I is given 



ical parameters and the Ginzburg temperature for a number of models, in particular, the 
pure solvophobic model, the pure Coulombic model as well as for the models including both 
the Coulomb interactions and the solvophobic interactions. In the present paper, all quan- 
tities, the critical parameters, the relevant coefficients of the effective Hamiltonian, and, 
accordingly, the Ginzburg temperature, for all the models considered have been calculated 
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in the one-loop approximation. This approximation leads to the mean-field values for the 
reduced critical parameters. It is worth noting that the charge-asymmetry does not manifest 
itself at this level of approximation. 

Having excluded the Coulomb interaction, we first calculate the reduced Ginzburg tem- 
perature for a simple fiuid model, the one-component system of hard spheres interacting via 
the square-well potential of the range A = 1.5. This well-known model is used to calibrate 
the Ginzburg temperature for the models having Coulomb interactions. Then, neglecting 
the short-range attraction we study the criticality in the RPM. Our calculations have shown 
that the reduced Ginzburg temperature of the Coulombic systems is much smaller (although 
it is not extremely small) than for simple fiuids. We have obtained the fiuctuation domi- 
nated region for the RPM to be much narrower (reduced by a factor of ~ 20) than for the 
HCSW fiuid. 

We have also studied the Ginzburg temperature depending on the interplay between the 
Coulombic and solvophobic interactions. Having introduced the ratio a that determines the 
strength of the Coulomb interaction with respect to the solvophobic interaction (SW poten- 
tial), we calculate both the reduced critical temperature and critical density, the Hamiltonian 
coefficients and the reduced Ginzburg temperature for a set of the a values ranging from 
0.5 to 100. For the model under consideration, the Ginzburg temperature shows a non- 
monotonous behaviour with the variation of a. In particular, Iq first decreases approaching 
its minimum in the region 10 < a < 25 and then starts to increase approaching the RPM 
value for a > 100. The present results generally agree with our previous results obtained 
within the framework of the different perturbative treatment. Moreover, they are in good 
quantitative agreement for a < 10. 

In conclusion, within the framework of the unified approach we have derived the Ising-like 
effective Hamiltonian which enables us to make reliable estimates of the reduced Ginzburg 
temperature for the RPM as well as for the model with the competing solvophobic and 
Coulomb interactions. For the model including both the short-range attractive interac- 
tions and the long-range Coulomb interactions, we have shown that the reduced Ginzburg 
temperature approaches the value obtained for the RPM when the strength of the Coulomb 
interactions becomes large enough. This suggests the key role of the Coulomb interactions in 
the crossover behaviour observed experimentally in ionic fluids. We believe that our results 
provide new insights into the nature of the non-classical region in the Coulomb-dominated 
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systems. 

VI. APPENDIX 

A. Expressions for coefficients ail" 



Taking into account the recurrence formulae for Tin (see (46) in Ref. 30|) we get the 
following expressions for coefficients a„ for n < A 

'J2,HSlfcj %Z 

M _ '^3,HS (2) l_ (3) _ (1-^) 

tin o 1 tin 



'S'a.HS 


03 ' 
'-'2,HS 

1 / 


o4 1 
'^2,HS \ 

_ 2(1 -^) 



^ " giz' ' ~ <ilz' 



2 



(0) _ ^ / C ^'^3,HS \ (2) 



4,HS - ^ I , "4 - -27' 

(3) _ 2{l-z) (4) _ 2{1- z + z-") 

^' ~ qlz^ ' ""^ " go'-2' ' 

where 5'„,hs = G'„_hs/(^)hS) Gn,ns is the Fourier transform of the ra-particle connected 

correlation function of a one-component hard-sphere system. 

In order to get the corrections Aas and Aa4 (Eqs. (l20l)-(l2Tl)) we also need coefficients a\ 

and Og presented by 

(2) _ _ 3! (2) _ 4! 
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